Influence of nanoconfinement on the nematic behavior of liquid crystals 
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We explore the nematic ordering of the rod-like liquid crystals 5CB and 6CB, embedded into 
parallel-aligned nanochannels in mesoporous silicon and silica membranes as a function of mean 
channel radius (4.7 < R < 8.3 nm), and thus geometrical confinement strength, by optical birefrin- 
gence measurements in the infrared region. The orientational order inside the nanochannels results 
in an excess birefringence, which is proportional to the nematic order parameter. It evolves contin- 
uously upon cooling with a precursor behavior, typical of a paranematic state at high temperatures. 
These observations are compared with the bulk behavior and analyzed within a phenomenological 
model. Such approach indicates that the strength of the nematic ordering fields a is beyond a 
critical threshold o c = 1/2, that separates discontinuous from continuous paranematic-to-nematic 
behavior. In agreement with the predictions of the phenomenological approach a linear dependency 
of a on the inverse channel radius is found and we can infer therefrom the critical channel radii, 
R c , separating continuous from discontinuous paranematic-to-isotropic behavior, for 5CB (12.1 nm) 
and 6CB (14.0 nm). Our analysis suggests that the tangential anchoring at the channel walls is of 
similar strength in mesoporous silicon and mesoporous silica membranes. A comparison with the 
bulk phase behavior reveals that the nematic order in nanoconfinement is significantly affected by 
channel wall roughness leading to a reduction of the effective nematic ordering. 

PACS numbers: 64.70.pm, 77. 84.Lf, 78.67.Rb 



INTRODUCTION 

The study of spatial confinement effects on liquid crys- 
tals (LCs) embedded in nanoporous and mesoporous ma- 
terials has been a subject of numerous experimental and 
theoretical investigations in the past[l-15|. It contin- 
ues to attract considerable interest in many areas of or- 
ganic electronics (l6j , not to mention the widespread use 
of such hybrid material systems in display technology. 
Memory effects resulting from frustration, as explored in 
recent computer simulations [ItJ for nematic LCs con- 
fined in bicontinuous porous structures, may be consid- 
ered as one prominent example where orientationally or- 
dered materials are expected to exhibit surprising func- 
tionalities because of topological confinement. However, 
implementation of such ideas in devices based on nanos- 
tructured materials requires a better understanding of ge- 
ometrical confinement and interface coupling that affect 
spontaneous nematic ordering. This problem has been 
emphasized in recent molecular dynamic simulations [l8| . 

The unusual behavior of confined LCs is often traced 
to pronounced quenched disorder effects which have been 
experimentally explored either by embedding of LCs to 
porous substrates [p-lll or by dispersion of small particles 
in LC fluids 0, ll"oj |. Strongly modified phase behavior 
occurs usually when the system is confined on submi- 
crometer length scales, where the characteristic size of 
pores or aerosil particles is just a few tens of molecular 
diameters. In this case the nematic ordering does not de- 
velop via a first-order phase transition, typical of the bulk 
state, but seems to be continuous with an orientational 



correlation length never increasing beyond the pore size 
Consequently, at temperatures even far above the 
bulk isotropic-to-nematic phase transition there exists a 
weak residual nematic ordering located preferably in the 
interface region (see Fig. 1(b)). Therefore it is more ap- 
propriate to term the corresponding phase as "parane- 
matic" rather than "isotropic" {f|. A large number of 
publications concern the nematogen alkylcyanobiphenyl 
LC family (nCB, n=5,7-9,12) embedded in a variety of 
porous solid substrates, including the random intercon- 
nected pore network of controlled pore glass matrices 
H| and Vycor glasses @ or the parallel-aligned cylindrical 
pores of aluminum oxide (Anopore) [3Lll9|, polycarbon- 
ate Nuclepore [2(J and porous silicon [J. The results ap- 
pear to be strongly dependent on the average pore size, 
types of LC investigated and porous matrix as well as 
the chemistry (surface-treatment) of the pore walls 0]. 
Indeed, these studies corroborate that the general behav- 
ior of nanoconfined LCs is concomitantly dominated by 
the anchoring force and the geometric confinement. On 
a phenomenological level they are described by a single 
quantity, the so-called nematic ordering field a. Relevant 
concepts have been developed by Kutnjak, Kralj, Lahaj- 
nar, and Zumer (KKLZ model) According to their 

considerations the nematic ordering field is proportional 
to the anchoring strength W and inversely proportional 
to the pore radius R. With increase of the ordering field 
the nematic ordering becomes more pronounced as it is 
sketched in Fig. 1(e). This effect is shown to be similar 
to that of an external magnetic field on a spin system 
which replaces a discrete phase transition by a contin- 
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FIG. 1: (color online). Orientational ordering of a nemato- 
gen LC inside a cylindric channel. Panels (a) and (b) sketch 
the nematic (T < Tj N ) and paranematic (T > Tj N ) order- 
ing, respectively, that develop in channels with smooth walls. 
Panels (c) and (d) depict the nematic ordering in broad and 
narrow channel with rough pore wall, respectively; dash-dot 
rectangles (red online color) mark the core region of the chan- 
nel filling, with relatively well developed nematic ordering 
in comparison to the considerably orientationally disordered 
near-interface region. Panel (e) shows the order parameter q e 
in the KKLZ model as a function of reduced temperature t for 
several effective surface fields a; traces 1-4,4' correspond to a 
case of smooth channel wall; channel wall roughness reduces 
the effective ordering as sketched by trace 4". 



uous evolution of the order parameter 0]. Experimen- 
tally, such a behavior has been originally demonstrated 
by NMR measurements on 5CB and 8CB confined to 
controlled porous glasses of different characteristic void 
sizes R @. However, due to the insufficient accuracy of 
the experimental data this analysis has been restricted 
to semi-quantitative evaluations only. Further progress 
has been achieved in recent high-resolution optical bire- 
fringence study on rodlike nematogens 7CB and 8CBfllT| 
confined in an array of parallel-aligned nanochannels of 
monolithic silica membranes. These matrices turn out 
to be particularly suitable for a very precise determi- 
nation of the order parameter behavior. However, the 



measurements in Ref. have been performed for one 
channel radius only (2i?=10 nm), thus a challenging issue 
remains a verification of the KKLZ model, particularly 
the i?-dependence of the ordering field a in the nanoscale 
region. 

In this paper we explore the nematic order- 
ing of the model nematogen LCs 5CB (4-Cyano- 
4'-pentylbiphcnyl) and 6CB (4-cyano-4'-hcxylbiphenyl) 
embedded in parallel-aligned nanochannels of size- 
controlled porous silicon (pSi) and porous silica (pSiC^) 
membranes studied by optical birefringence measure- 
ments. The experimental results obtained are compared 
with the bulk behavior of these LCs and analyzed within 
the KKLZ approach. The choice of these low-molecular 
LC compounds is motivated for the two following rea- 
sons. First, in the bulk state both thermotropic liquids 
are characterized by quite simple phase diagrams. Upon 
cooling they exhibit a single phase transition before solid- 
ification, i.e. a transition from the isotropic (I) to the ne- 
matic (N) phase at Tin ~ 308.5 K (5CB) and T IN «302.4 



K (6CB) [2l|, |22j, respectively. Second, the strength of 
the first-order character of the isotropic-to-nematic tran- 
sition practically is not affected by the confinement as 
has been ascertained in recent theoretical studies 23, 24 1 . 



The latter essentially simplifies the analysis based on the 
KKLZ model. 



EXPERIMENTAL 

Nematogen LCs 6CB and 5CB have been purchased 
from Merck. The pSi membranes have been prepared by 
electrochemical anodic etching of highly p-doped, (100) 
oriented silicon wafers. To obtain membranes of pSiC>2, 
the free standing pSi membranes have been subjected to 
further thermal oxidation for 12 h at T=800 °C under 
standard atmosphere. By applying different etching con- 
ditions we have obtained membranes with average pore 
radius 4.7±0.3 nm (pSi, porosity P=60%), 6.3±0.3 nm 
(pSi0 2 , P=50%) and 8.3±0.3 nm (pSi, P=67%) as deter- 
mined by recording of volumetric N2-sorption isotherms 
at T=77 K. The mesoporous membranes were completely 
filled by capillary action (spontaneous imbibition) of LCs 
in the isotropic phase 25(. For the bulk measurements 
the sample cells, made of parallel glass plates (d=30 



/im), have been filled by the LCs providing homogeneous 
(5CB) or homeotropical (6CB) alignments. 

Note that pSi- membranes of about 300 /xm thickness, 
as used in our studies, are not suitable for optical polar- 
ization measurements in the visible spectral region due 
to a strong fundamental light absorption at the wave- 
lengths below ^1000 nm. For this reason the optical 
birefringence measurements were performed at A =1342 
nm in the infrared region, for which both pSi and pSiC>2 
membranes, the glass sample cell employed for the bulk 
measurements and both host LC compounds, are opti- 
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cally well transparent. 

The samples were tilted out with respect to their long 
channel axis (optical axis) by an angle of 40-45 deg. The 
optical polarization setup employs a photoelastic modu- 
lator and a dual lock-in detection system, as described in 
Ref . [n[ , providing an accuracy of the optical retardation 
measurements better than 5T0 -3 deg. 



RESULTS AND DISCUSSION 

In Fig. 2 we show the temperature dependences of the 
optical birefringence measured for the mesoporous mem- 
branes with different pore radii filled by LCs 5CB and 
6CB, respectively. The nematic ordering does not oc- 
cur via discontinuous, first-order phase transitions, as 
observed in the bulk state (see Fig. 3), but evidently 
evolves continuously in the vicinity of the temperature 
corresponding to the isotropic-to-nematic bulk transition 
Tin. 

Consequently, at temperatures above Tin there is a 
specific precursor behavior indicating a residual nematic 
ordering at high temperature, typical of a paranematic 
state. The continuous character of the paranematic-to- 
isotropic transition can be described by a sufficiently 
large nematic ordering field a > a c — 0.5 (due to the 
small channel diameters) , if one refers to the KKLZ ap- 
proach. By the same token a comparison of our mea- 
surements with the predictions of this model allows us to 
determine a for each channel radius R as outlined below. 

The untreated surface of silica and silicon pores en- 
forces planar anchoring without a preferred lateral direc- 
tion (within the plane) [3, [26j], whereas the elongated 
geometry of cylindrical pores yields additionally a pre- 
ferred orientation of molecules along the pores axes. In 
accordance to the KKLZ model it is convenient to char- 
acterize the relevant ordering by the scaled order param- 
eter q = Q/Qo, where Q = ^(3 cos 2 9 — 1) stands for the 
nematic (orientational) order parameter in its classical 
definition, Qq = Q{Tin) and is the angle between the 
long molecular axis and the director. The dimensionless 
free energy density of a nematic LC in confined geometry 
reads then: 



f = tq 2 - 2q 3 + q i -qa + nq 2 



(1) 



where t = {T - T*) / {Tin - T*) is the dimension- 
less reduced temperature, T* is the effective tempera- 
ture and the K-term completes the free energy expansion 
0,111 by a contribution originating from quenched dis- 
order effects. Minimization of Eq.([T]) with respect to q 
gives its equilibrium magnitude q e as it is sketched in 
FigHJe). The isotropic (paranematic)-to-nematic tran- 
sition occurs at t n = 1 + a — K and is discontinuous as 
long as fj < <t c = 0.5. Above this value a continuous 
transition is realized. 



The orientational order inside the pores results in an 
excess birefringence, 6 (An) oc q e , which appears on the 
background of the geometric birefringence, An 9 , be- 
ing a characteristic feature of membranes with paral- 
lel nanochannels (27j . The measured birefringence rep- 
resents a superposition of these two contributions, i.e. 
An = 5 (An) + An g . The geometric birefringence is quite 
small for the entirely filled pSi02-membrane due to simi- 
lar refractive indices of host silica (n re 1.46) and isotropic 
(or average) refractive indices of the guest 5CB or 6CB 
LCs (n rel.58 12110 . For this reason it has been ignored in 
earlier studies [11] . However, the geometric birefringence 
is considerable in the case of entirely filled pSi membranes 
because of the large refractive index of the silicon host 
(n re3.4 [111). Thus 6 (An) and An g appear to be com- 
parable, as it is evident from Fig.[2Jc,f). Moreover, An g 
exhibits a smooth temperature behavior (referred to as 
" thermooptic" changes in the following) due to the vary- 
ing refractive indices of the host Si (Si02) matrix and 
the guest LC as well as changes of the porosity caused 
by thermal expansion/contraction of the porous matrix. 
Accordingly, this "thermooptic" contribution, approxi- 
mated by a linear dependence, has been included in the 
fitting procedure, see grey dash-dot line in Fig. [2] The 
fitting procedure is based on an iterative, numerical min- 
imization of Eq.(l) to extract a and k values that pro- 
vide smallest least-square deviations between the calcu- 
lated curve and the experimental points. Similarly to 
Ref.0 we first calibrated the temperature scale of the 
KKLZ model by analyzing the bulk birefringence in the 
region of the isotropic-to-nematic transition. This anal- 
ysis gives Tin ~ T* approximately equal 3.4 K (5CB) 
and 3.2 K (6CB), resp., i.e. values close to the magni- 
tude of 3.5 K obtained for 8CB 0. In fact the difference 
Tin—T*, which defines the proximity of the system to the 
tricritical point, may vary in general under confinement 
conditions, as has been recently emphasized by Kralj et 
al 23J, |24| . This effect appears to be dominant if surface 
wetting interactions are weak, which is observed in sam- 
ples containing relatively highly flexible LCs molecules, 
as e.g. 12CB. However, it is expected to be weak in 5CB 
and 6CB, where surface wetting strengths play impor- 
tant roles. This particularly justifies a choice of these 
LCs for the present experiments. Because of this rea- 
son the temperature conversion factor, defined by the 
difference Tin ~T* , has been fixed in all consequent fit- 
ting procedures to values specified above. Dashed red 
curves in these figures represent the best fits of measured 
An(T)-dependences based on the free energy expansion 
([T]). Including higher order terms (up to the tenth or- 
der) into Eq. ([l} better reproduces the saturation of the 
nematic order parameter at low temperatures (solid red 
lines in Fig. [2]), but has no influence on the fit quality 
in the transition region or the extracted a and k val- 
ues. Insignificant deviations between the experimental 
points and fitting curves may be recognized also directly 
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FIG. 2: (color online). Measured optical birefringence An (A =1342 nm) of nematogen 5CB (left panels) and 6CB (right 
panels) embedded in parallel-aligned nanochannels of several porous membranes: (a, d) pSi (P=68%, R = 8.3 nm); (b, e) 
pSi02 (P=50%, R — 6.3 nm); (c,f) pSi (P—60%, R = 4.7 nm). The broken line (red online color) represents the best fit 
obtained within the KKLZ- model based on a free energy expansion {TJ. Completion of the free energy by higher order terms 
(up to tenth order) describes the saturation of the order parameter at lower temperatures more accurately, see solid line (red 
online color). The base line corresponding to the thermooptic contribution (dash-dot line, grey online color) is obtained within 
the same fitting procedure. The values of the effective nematic ordering field a and quenched disorder factor k are quoted for 
each LC/porous matrix combination. 



in the transition region. In most cases the fits exhibit 
here slightly steeper changes than the experimental data. 
The reason for such discrepancies could originate in pore 
wall roughness as well as the pore diameter distribution, 
which arc both not considered in the KKLZ model and 
which should lead to a more smeared behavior in the 
transition region. 

In Fig. [3] we show a vs It is for both confined LCs 
investigated practically linear. Following the considera- 
tions in Refs. 0, H| the effective nematic ordering field a 
is expressed as: 



RkQo 



(2) 



where £* = £(Tjjv) is the correlation length at the tem- 
perature of the isotropic-to-nematic transition, W n \ is 
the anchoring strength and k is the nematic elastic con- 



stant. Note that we neglect here the first term of Eq. 
(7c) of Ref.[|| which is oc R~ 2 . In the case of tangential 
ordering it is expected to be small. For this reason the 
ct(.R _ ^-dependence in Fig. [3] is dominated for each LC 
by the linear contribution in accordance with Eq. @ . In- 
terestingly, for both LCs the data point <r(R =6.3 nm) 
determined for the pSiC>2 membrane lies perfectly on the 
linear fit which passes through two other data points, 
a(R =4.7 nm) and a(R =8.3 nm) that correspond to pS\ 
membranes. This implies an identical anchoring strength 
for the silicon and silica membrane, which appears sur- 
prising at first glance. We suggest that the native oxide 
(SiC^) layer usually forming on silicon surfaces renders 
the inner surfaces of the silicon membrane silica-like and 
thus the fluid/ wall interaction, which determines the an- 
choring strength, is identical in porous silicon and silica. 
Slightly different slopes of the tr(i? _1 )-dependences, as 
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FIG. 3: (color online). Nematic ordering field a vs inverse 
mean channel radius R -1 . Square symbols correspond to a 
values as extracted within the fitting analysis of the An(T)- 
dependences shown in Fig. [2] Solid lines are the best linear 
fits of the a data points. Their intersections with a horizontal 
line, a = a c = 1/2 (see solid circles) yield coordinates of the 
critical points a c (R c = 12.1 nm) (5CB) and a c (R c = 14.0 
nm) (6CB) as discussed in the text. 



one can see in Fig. [3] for nanoconfined 5CB and 6CB 
LCs, originate from the difference in material constants 
characterizing these LCs and, likewise, the difference in 
anchoring forces. The slope of the cr(i? _1 )-dependence 
is denned by the factor g = 2£* 2 VF„i/(fcQ ). The fit 
of the data points in Fig. [3] yields g-factor of 6.05 nm 
(5CB) and 7.03 nm (6GB), respectively. With the typ- 
ical correlation length £* ~10-15 nm, clastic constant 
k =7.7-10- 12 N (5CB), k =8.0-10' 12 N (6CB) [U and 
Qo ~l/2 one yields an anchoring strength W n \ in the 
range 50-100 p;J/m 2 . These values arc of the same order 
of magnitude as obtained for 8CB in random pore net- 
works of controlled pore glass Q, i.e. the interaction of 
5CB or 6CB molecules with SiC>2 surface corresponds to 
a regime of weak anchoring. 

Interpolation of the a(R~ 1 ) dependences to er c =0.5 
(see Figj3|) allows one to determine the critical pore radii 
R c 12.1 nm for 5CB and 14.0 nm for 6CB, respectively, 
see solid circles at the intersection of relevant lines. These 
values arc close to the ones obtained for other nCB LCs, 
e.g. 7CB (R c =8 nm) or 8CB (i? c =11.5 nm) embedded 
into pSi02 membranes (llj . 

By subtracting the thermooptic contribution, An g {T), 
from the measured birefringence An(T) shown in Fig. [2] 
one obtains the excess birefringence, 8(An) oc q e , origi- 
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FIG. 4: (color online) . Normalized excess birefringence An nex 
vs temperature T (blue, green and red colors) of nematogen 
5CB [panel (a)] and 6CB [panel (b)] embedded into meso- 
porous membranes, as calculated from the data presented in 
Fig. [2] The traces An nex (T) for pSi0 2 (6.3 nm, blue color) 
and pSi (8.3 nm, green color) are shifted vertically by 0.015 
and 0.03, respectively, in order to avoid their overlap; hor- 
izontal broken lines with relevant labels marks 0-levels for 
these traces. Bulk birefringence (black circles) is given for 
comparison. 



nating exclusively from the orientational ordering of the 
confined molecules. The excess quantity 5(An) roughly 
scales with the porosity of the membrane, P, thus a com- 
parison of the orientational ordering in samples of differ- 
ent porosity or with the bulk state requires its normal- 
ization by P. In a more accurate characterization, which 
is based on the Bruggeman approximation, the normal- 
ized excess birefringence is given by An nex = (5(An)/P, 
where £ is a correction factor that depends on composite 
parameters, i.e. the actual pore geometry, the porosity 
and the difference between the refractive indices of the 
host material and the guest LC. If the latter is small, as 
e.g. for the composites pSi02+nCB (n=5-10), the factor 
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FIG. 5: (color online). An nex / An ( - max ^ vs _R _1 as determined 
for both LCs at T — Tf N — 10 K using the normalized excess 
birefringence, An nex (T), and bulk birefringence, An(T), data 
in Fig. H 



C deviates only very slightly from 1.0 (£ =0.993-0.994) 
and the corresponding correction can be ignored. This, 
however, is not the case for the composites with consid- 
erably different refractive indices of the host and guest 
components as, e.g., pSi+nCB for which ( equals 0.943 
(P=60%) and 0.909 (P=68%) as provided by estima- 
tions based on the Bruggeman approach for membranes 
with strongly elongated parallel channels; for details re- 
garding the methodology of such evaluations see Ref. [27| ■ 
The normalized excess birefringence An nex of nanocon- 
fined LCs, determined in such a way, is shown in Fig. 0] 
and is compared with the bulk behavior (black color). 

In general, a decreasing channel radius leads to an in- 
crease of the effective nematic ordering field a. Accord- 
ingly, enhancing of the nematic ordering is expected in 
the entire temperature region as shown in Fig. [He); see 
the sequence of solid green, blue and red lines correspond- 
ing to a rising sigma value (k~0; 02 < 03 < 04, respec- 
tively). On the other hand, quenched disorder effects, 
represented by the K-term in the free energy (TjQ) lead to 
a lowering of the paranematic-to-nematic transition as it 
is illustrated by the dashed red line in Fig. [1] On the 
reduced temperature scale it leads just to a temperature 
shift of the whole curve by At = —k. Thus depending 
on the relative magnitude of n and a the orientational 
ordering at a certain fixed temperature below Tin may 
even decrease while the effective nematic ordering field 
rises. 

Distortions caused by surface roughness of the pore 
walls, as sketched in Fig.[TJc)(d), lead to further decrease 
in the effective (averaged) nematic ordering. Accord- 
ingly, the dashed curve 4' in Fig. HJc) in fact corresponds 
to the theoretically achievable maximum value qi m 
for, so to say, ideally smooth channel walls [Fig. [IJb)], 
whereas the dash-dot curve 4" accounts for disordering 



due to channel wall roughness [Fig. [He), (d)] leading thus 
to its reduced magnitude q*. By 

An (max) we 

denote the 

birefringence relevant to qe , it can be evaluated for 
each pore geometry and temperature using the measured 
bulk birefringence (Fig. @] (a)(b)) and the extracted fit 
parameters, a and k. On the other hand, An nex is rel- 
evant to q*. In such a consideration the ratio q*/qi max ^ 
characterizes the degree of confined nematic order and is 
equal to the ratio An nex / Arv- max \ In Fig. [5] we present 
An nex / ArS max ^ vs as determined for both LCs at 
T = Tjn — 10 K. The effective orientational order evi- 
dently rises with increasing pore diameter. Such a trend 
can be rationalized by a consideration of the contribu- 
tions that originate from different regions of a pore fill- 
ing where relatively well developed nematic order in the 
core region competes with orientational disorder of few 
molecular layers next to the pore walls, as sketched in 
Fig. [H(c) an d (d). Assuming that the nematic order is 
fully developed in the core region and is negligible in the 
near-interface region the degree of nematic order inside 
the pore q*/qi max ^ = (1 — Iq/R) 2 , where lo is the ef- 
fective thickness of the disordered near-interface layer. 
Here Iq is of the order of few molecular layers satisfy- 
ing the inequality, < Ro- For large pore radii the 
effective nematic order is favored by the ordered core 
component {cfe/qe \r^,oo — 1 ), whereas at small ones 
the disordered near-interface layers provide a dominant 
contribution {ql/qi max ^ \r->i = ). The dashed black 
line in Fig. [5] represents the best fit of the experimental 
data points by the (1 — l / R) 2 -hmction which yields an 
effective thickness of the disordered near-interface layer 
of Iq =1.35 nm, corresponding to approx. 2 to 3 inter- 
molecular spacings ao (for 5CB a = 5.0 A 30]). One 
should emphasize that represents here the effective 
model parameter rather than the thickness of the real 
near-interface layer. In reality the near-interface region 
exhibits a certain nematic ordering. Thus, presumably, 
our case corresponds to much weaker distortions of the 
nematic order that start in the layer next to the pore 
wall and radially decays on distances comparable with 
the orientational correlation length, i.e. at least over dis- 
tances several times longer than Iq. Such a radial gradi- 
ent in the orientational order has been inferred in Monte 
Carlo simulations on rod-like liquid crystals confined in 
nanochannels 31, 331 ■ 



4. CONCLUSION 

We have reported an optical birefringence study with 
infrared Laser light on the nematic ordering of the 
nematogen LCs 5CB and 6CB, embedded in parallel- 
aligned nanochannels of channel-radius-controllcd porous 
silicon and silica membranes (4.7 < R < 8.3 nm). These 
results are compared with the bulk behavior and ana- 
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lyzed within a phcnomcnological approach. 

The molecular oricntational order inside the channels 
gives rise to an excess birefringence, which is proportional 
to the nematic order parameter and appears on top of the 
background of the geometric birefringence typical of a 
homogeneous matrix permeated by parallel-aligned nano 
channels. By applying an appropriate fitting procedure 
we were able to separate these distinct birefringence con- 
tributions and thus to determine the magnitude of the 
effective nematic ordering field a originating from the 
geometric confinement. 

The nematic ordering in the cylindrical nanochanncls 
evolves continuously with characteristic precursor behav- 
ior typical of the paranematic state. The strength of the 
ordering fields, imposed by the silicon or silica walls, ap- 
pears in all cases investigated beyond a critical threshold 
a c = 1/2 , that separates discontinuous from continu- 
ous paranematic-to-nematic behavior. An interpolation 
of the linear cr(i? _1 ) -dependence provides values for the 
critical pore radius, R c = 12.1 nm (5CB) and R c = 14.0 
nm (6CB). Those values are close to the ones obtained 
recently for 7CB and 8CB embedded int porous silica 
membranes. 

Our analysis evidently suggests that the tangential an- 
choring strength acting at the channel walls is the same 
in porous silicon and silica membranes, presumably be- 
cause of the formation of a native oxide layer at the inner 
surfaces of the porous silicon membrane. A comparison 
with the bulk behavior suggests that the nematic order 
of 5CB and 6CB in nanoconfinement is significantly af- 
fected by the roughness of the channel walls. The impact 
appears to be more pronounced for narrow pores (R < 
5 nm), where the reduction in nematic ordering exceeds 
even 50%. 
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